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| The spontaneous supersymmetry-breaking that takes place in certain spin-glass models signals a 

particular fragility in the structure of metastable states of such systems. This fragility is due to 
the presence of at least one marginal mode in the Hessian of the free energy, that makes the states 
highly susceptible under external perturbations. The cavity method is a technique that recursively 
describes the property of a system with TV + 1 spins in terms of those of a system with TV spins. To 
do so, the cavity method assumes a certain degree of stability when adding a new spin to the system, 
i.e. it assumes that for a generic choice of the parameters there is an one-to-one correspondence 
between the metastable states of the system with TV spins and the metastable states of the system 
with TV + 1 spins. In systems where the supersymmetry is broken such a correspondence does 
not exist, and an alternative formulation of the cavity method must be devised. We introduce a 
generalized cavity approach that takes care of this problem and we apply it to the computation of 
the probability distribution of the local magnetizations in the Sherrington-Kirkpatrick model. Our 
findings agree with the correct supersymmetry-breaking result. 
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In physics we often meet two quite different situations. On the one hand, we tend to assume (or hope) that small 
perturbations have small effects. This is at the basis of all attempts to use any sort of perturbation theory. On the 
other hand, we know that many interesting phenomena take place when such an assumption is in fact violated. At 
criticality, whenever the susceptibility of a system is infinite, its response is anomalous. In this case small perturbations 
may indeed have big effects. This state of affairs is often accompanied by the presence of zero, or marginal, modes 
among the second derivatives of the action, may these be the masses, or the inverse susceptibilities. When this 
happens we must be extremely careful in applying perturbation theory. 

A particularly interesting realization of this scenario occurs when criticality is induced by spontaneous symmetry 
breaking. In that case, symmetry breaking, marginality, and break-down of standard perturbative techniques, give 
rise to some interesting phenomena. In this work we discuss how the spontaneous breaking of a supersymmetry is 
. connected to the presence of marginal modes and anomalous response in certain spin-glass models. The main focus of 
our study is the structure of metastable states in such models, and in particular the extreme fragility of this structure 
under small perturbations. A significant consequence of this feature is the break-down of the standard cavity method 
0,0, whose basic assumption is that by adding one extra degree of freedom to a large systems, its physical properties 
(and in particular some properties of the structure of states) do not change dramatically. This assumption is no longer 
valid when the supersymmetry is broken. Our main result is to provide a generalization of the cavity method that 
. works also in supersymmetry-breaking systems. 

Metastable states in mean-field spin-glasses can be identified with the local minima of a mean-field free energy F 
(also known as the TAP free energy 3]), that is a function of the local magnetizations m t of the system. Generally 
speaking, in order to compute the number of local minima of a function F(m) , one can introduce an effective action 
that is invariant under a generalized form of the Becchi-Rouet-Stora-Tyutin (BRST) supersymmetry [J, |fj, |6|, |7], |8| • 
The Ward identities generated by this symmetry have a rather clear physical meaning. The most relevant one reads 
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where Aij(m) = [didjF(xn)] 1 is the inverse of the second derivative, i.e. of the Hessian. The brackets (•) indicate 
a sum over all metastable states of the system. The meaning of equation JQ) is straightforward: it expresses the 
natural relationship between the susceptibility and the curvature of the minima of the free energy. Therefore, the 
supersymmetry seems to encode a very robust physical feature of metastable states, namely the static fluctuation- 
dissipation theorem. This notwithstanding, it is now believed that in certain spin-glass models the supersymmetry 
is in fact spontaneously broken in the low-temperature phase, and that relation Q is thus violated |3, lift 1 1 1| . Such 
systems include the Sherrington-Kirkpatrick (SK) model 0], that we discuss here. 

The reason why the static fluctuation-dissipation theorem is violated, and thus the supersymmetry is broken, lies in 
the peculiar structure of metastable states of such systems. Recent studies show that at low temperatures 
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all stationary points of the mean-field free energy are organized into minimum-saddle pairs. The minimum and the 
saddle are connected along a mode that is softer the larger the system size N. Moreover, the free energy difference of 
the paired stationary points decreases with increasing A. In other words, supersymmetry-breaking metastable states 
are marginal in the thermodynamic limit, having at least one flat zero mode. In this situation it is clear that even an 
infinitesimal external field may destabilize some states, making them disappear. On the other hand, virtual states, i.e. 
inflection points of the free energy with a very small second derivative, may be stabilized by the field, giving rise to 
pairs of new states. In such a situation we must reconsider the validity of equation Q). At the l.h.s. we differentiate 
with respect to an external field a sum over all metastable states. The problem is that, due to marginality, some 
elements in this sum may disappear or appear as the field goes to zero. Therefore, even though the static fluctuation- 
dissipation relation holds for each individual state, when we sum over all states, an anomalous contribution arises due 
to the instability of the whole structure with respect to the field. Relation is thus violated. 

Supersymmetry breaking is thus the mathematical expression of a great instability in the structure of metastable 
states. Adding a new spin to the system, the basic building block of the cavity method, will therefore not be harmless. 
Our aim is to find a more general formulation of the method, valid also in supersymmetry-breaking systems. Note that 
a similar task has been independently pursued in 0] , although introducing an explicit violation of the supersymmetry. 
Our starting point is the set of equations satisfied by the local magnetizations in a system with A degrees of freedom. 
These are just the equations of stationarity of the free energy, 
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£T J (m) = i=l,...,N, (2) 



dm. 

where m = mi, . . . , mjv. For the SK model we have Q, 

N 

e\ N] = tanh" 1 m 4 + 1 (1 - q) m, - f3 ^ J ljmj . (3) 
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The self-overlap is q = (l/A^J^Vmf, an( A the random couplings Jy are drawn from a Gaussian distribution of mean 
zero and variance 1/A. Let us introduce the number dAf( N >(m) of solutions of eqs.(j2Jl with magnetizations between 
m and m + dm. The density of solutions pW(m) is then defined as, 

dyVW(m) =pW(m) dm . (4) 

The basic point of the cavity method 0,0] is to write a recursive relation expressing the density p( N+1 ) of the system 
with A + 1 spin, as a function of the density of the system with A spins. The fundamental idea behind this 
procedure is that adding one single spin to a system with a large number of degrees of freedom, is a small perturbation 
and consequently the properties of the new A + 1 system will not be much different from the original A system. More 
specifically, a natural assumption is that the structure of metastable states of the system does not change when we add 
the new spin: the solutions of @ are slightly modified, but there is still a one-to-one correspondence between solutions 
of the A and A + 1 systems. This hypothesis seems reasonable, and it has been at the basis of all standard cavity 
investigations. However, we will show that it is in fact equivalent to assume that the supersymmetry is unbroken. 

Let us add a new spin at site 0, and call (m, mo) the global magnetization vector of the A + 1 system. The new 
set of A + 1 equations can be split as, 

El N+1 \m,m ) = i = l,...,N, (5) 

E (m,m ) = . (6) 

We can examine this set of equations by first solving the first A equations © at fixed mo, and then plugging m(mo) 
into ©. The assumption of stability of states can be expressed by saying that there is a one-to-one mapping between 
solutions m^ 7V+1 - ) of eqs.(JHJ) and solutions of eqs.©: if is a solution of the A system, the effect of adding the 
spin is just to slightly modify the old magnetization to a new value m^ A ' +1 ^ = f mo (m^). From this key hypothesis, 
it follows that for a given mo the number of solutions with magnetization m( JV+1 ) of the first A equations of the 
A + 1 system, is simply equal to the number of solutions of the original A system that have magnetization m^), 

dA A(^+i) (m (Jv+i) ) = dAfW( m W) , (7) 

with m^ A ' +1 ^ = { mo (m( N )). In order to obtain the density of solutions of the new A + 1 system we must also impose 
the extra equation © for mo. This relation can be cast in a simple physical form, namely mo = tanh (PY]j Jojmi ): 



3 



the spin at zero aligns to the local field produced on site by the original N system. The density of solutions for the 
new system therefore becomes, (sums over repeated indices are understood), 

p (N+l) ^ m (W+l) j mo ) dm ("+D dmo = p (N) ^ m (N)^ dm (N) g ^ anh -l mQ _ pj 0jm <. N ^ _J_ dmQ , ( 8 ) 

The first part of this equation is equivalent to relation (JJJ , while the second part makes sure that the equation for 
mo is satisfied (the last factor is the Jacobian of the (5-function) . The meaning of this equation is clear: the density 
of solutions remains the same, provided that we satisfy the equation for the new site 0, and that we change a little 
the other N magnetizations. The map between the old and the new magnetizations can be found by noting that the 
new magnetization mo exerts a field Jiom$ on the old sites, and thus m\ N ~^ = m\ N ^ + \ Jw m 0i where x — — <z) 
is the susceptibility. 

Equation (JHJ) is a recursive relation for the number of states. In order to get a self-consistency equation, we have 
to work with the probability, i.e. the density divided by the total number of states, p(m) /Af, that has a well defined 
limit when N — > oo. In the annealed approximation we assume that the probability is self-averaging and that it can 
therefore be averaged over the disorder |l4j. The averaged probability density factorizes, i.e. p(m)/M — \\ i p{rrii), 
and thus Eq. JSJ becomes a self-consistency equation for the single site probability p(mj): we divide both sides by 
the global number of states M^ N+1 \ average over the disorder, and integrate over the first N magnetizations. In the 
thermodynamic limit p( N+1 ~> = p^ N \ and we finally get, 

p(mo) = / JJdP(Jio) drrn p{rrii) K, ss (m,m Q ) , (9) 
J i 

where, 

/C ss (m, mo) = S (tanh -1 mo — fiJojirij) /(l — tuq) . (10) 

The factor Af( N ^ /A/"*- Ar+1 * ) on the r.h.s. has been reabsorbed into the normalization constant of the probability p. This 
equation is not difficult to solve, and we find, 



p(mo) = FTffT exp 
V 27r P 1 



(tanh 1 mo) 2 
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(11) 



where q = (m 2 ), and where (X) = J dm p(m) X (m) . Equation (jl 1|> coincides with the probability distribution found 
in Q , that is the super symmetric distribution. The assumption of stability of metastable states is thus equivalent to 
assuming unbroken supersymmetry. However, the BRST supersymmetry is in fact broken in the SK model [jl Il0|. 
and thus the correct distribution is not given by l|ll|) . To find the supersymmetry-breaking distribution within the 
cavity approach, we must give up the assumption of stability of states, and find a more general formulation. The first 
step is to recognize that eqs.© for the N + 1 system can be formally rewritten as, 

£ 4 (JV+1) (m, moW) = E\ N \m\(i') - fc 4 (m,m ) = (12) 

with, 

ki(m,m ) = (3Jiom - — (1 - m 2 ,) m 4 . (13) 

The variance of Jij in the N + 1 system is smaller than in the N system. We have then rescaled the inverse temperature 
from P to j3' — (3[1 — \/{2N)] in such a way that the disorder appearing in e\ N ^ has the correct scaling for a system of 
size N. Equation (|12H suggests the path we have to follow: the function fc, may be seen as a local field acting on site i 
of the N system, so that, in a way, finding solutions of the N + 1 system is like finding solutions of the N system with 
a field. If this can be done, we may hope to write a self-consistency equation not simply for the probability density of 
solutions with given magnetization, but for a more complicated object, i.e. the probability density of solutions with 
given magnetization and field. The equations for the metastable states of an N system with a field are, 

f W^. (14) 
We define the new density of solutions p( N \m\h) as, 

aW {N) {m\h) = pW(m|h) dm , (15) 
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where dAf( N ^ (m|h) is the number of solutions of 114(1 . with given magnetization m and external field h. The idea is 
to write a recursive relation for p( N \m\h) rather than for p^ N ^(m). From (|12fl we have that the equations for the 
N + 1 system in an external field (h, h ) are, 

£f°(m) = ki(m,m ) + hi (16) 
E (m,m ) = h . (17) 

Let us focus now on l|16|) . at fixed value of too- The number of solutions of this set of equations seems formally equal 
to the number of solutions of a system with external field k + h. The only problem is that k is not an external field, 
since it depends on m, and thus d/V(m|h + k(m)) ^ p(m|h + k(m)) dm. Fortunately, it can be simply proved a 
similar relation, that reads, 

dA/" (Ar) (m|h + k(m)) = p {N) (m\h + k(m)) w(m) dm , (18) 
where the re- weighting factor is given by, 

d m [E^ N \m)] 

We can apply this general formula to equation i|16f) . and remember that (as in the previous case) to obtain the full 
solutions density for the N + 1 system we must also satisfy the extra equation for mo. This time, however, there is no 
mapping between solutions of the N and N + 1 systems, and thus the field experienced by spin must be expressed 
as a function of the new magnetizations Q , 

E = tanh -1 rn - j3 Joj m j + /? 2 (1 — l) m o , (20) 

j 

and we have to remember that m depends on too- In conclusion, we have, 

p (A,+1) (m,TO |h, h ) = p (Ar) (m|h + k(m, too)) w(m, to ) x 

x 5(E (m,m )-h ) — ^ . (21) 

CtTOo 

Equation (|21|l is a general recursive relation for the density of states, that is always valid, whether the structure of 
states is stable or not. The supersymmetric result can be recovered if we assume that there is a mapping between 
solutions m of the N + 1 system and solutions f~^(m) of the system in absence of the site. In this case we have 
(for any fixed N), 

p(m|h + k(m,TO )) w(m,m ) dm = p(f mo 1 (m)|h) df~]. (22) 

In other words, in the supersymmetric case the effective field k can be reabsorbed via a change of variables, giving 
back the supersymmetric expression (|5J). 

When we pass from single instances to probabilities, we can average equation H21|l over the disorder and once again 
write self-consistency equation for the single-site probability p(mi\hi). After some algebra we get, 

P/3(m \h ) = / ]^[d-P(J 4 o) drrii p^{mi\hi +ki(i7ii,m )) /C(m, m ) (23) 
•* i 

with, 

/C(m,TO ) = 5 [tanlT 1 to - /3J jTOj + (3 2 {l - q)m ~ ho] exp[/3 2 (l - m 2 Q ){\ - q)} - — — ^ , (24) 

1 — to 

and where we recall that, (3 1 = 0[1 — 1/(2N)]. The exponential in l|24() is the factor lu, while the last term is the 
Jacobian of the delta. Let us compare equation © to 123fl : the key point is that in the general case there is no closed 
equation for p(rrii), however complicated. On the other hand, if we start from a system with external magnetic field, 
the effect of the new spin can be counterbalanced by tuning the local fields in the new system, in such a way to leave 
stable the structure of states. For this reason it is possible to write a self-consistent equation only for the probability 
p(m,i\hi), but not for p(rrii). This is the core of our approach. 
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In order to solve equation l|23|) we adopt an ansatz for the asymptotic form of p(m\h) (we drop now the subscript 
0). For small external fields, we expect, 

p(m\h) = p(m) exp [A(m)h 2 + B{m)h] , (25) 

where terms of higher order in h have been discarded, and A(m) and B(m) are functions to be determined self- 
consistently. By setting hi = and h ~ ki ~ l/ViV in equation l|23J) . and by using an integral representation for the 
delta function, we get to leading order in N: 



p(m\h) = C exp 1 —z- [tanh 1 m — Am] 2 + Am 



CX P I — h 2 + 75T - [tanh 1 m — Ami /z I 
[ 2/3 2 g /3 2 g J 



(1 -m) 2 

(26) 



withg = (m 2 ), A = (A(m)>, A = (3 2 {mB(m)) - (3 2 (\-q) and A = {(3 2 /2)[(B(m) 2 ) + A] + A. We note that the shift in 
temperature only gives a constant term at the leading order: together with all other constant contributions it has been 
reabsorbed in the normalization constant C. By comparing (|25(l with (|26|) we immediately have, A(m) = — l/(2(3 2 q), 
B(m) = [tanh -1 m — Am]/(/3 2 g). The parameters A and A therefore satisfy the equations: 

1 P 2 

A = — (m tanh - m) (1 — q) 

2q x ' 2 K H ' 

A = A ~ Y q + W ([tanh_1 m " Am]2) ' (27) 

The expression of pirn) we have found (first line in equation l|26[0. and the equations for the parameters A and A 
are precisely those of the supersymmetry-breaking solution found long ago in |l5j . The supersymmetric expression is 
recovered for A = A = Q . 

In this work we have shown how the cavity method can be generalized to situations where the structure of states is 
unstable under external perturbations, i.e. where the supersymmetry is broken. In this case, the simple distribution 
p(m) is not stable when a new spin is added to the system. However, the perturbation caused by the new spin is 
analogous to an external field. Therefore, if we consider the system in presence of a field, we can balance the effect 
of the new spin by tuning the field. This means that the distribution p(m\h) is stable when the new spin is added, 
and a self-consistent relation can be written for it. It is important to stress that this equation can be written for any 
system, once the form of the effective field fc, is known. In particular, applications of this method to diluted systems 
are under study. We thank F. Ricci Tersenghi for many interesting discussions. 
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